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Abstract. This paper investigates the origin of the observed large variety in dust-to-gas ratio, T), among blue 
compact dwarf galaxies (BCDs). By applying our chemical evolution model, we find that the dust destruction 
can largely suppress the dust-to-gas ratio when the metallicity of a BCD reaches 12 + log(0/H) ~ 8, i.e., a 
typical metallicity level of BCDs. We also show that dust-to-gas ratio is largely varied owing to the change of dust 
destruction efficiency that has two effects: (i) a significant contribution of Type la supernovae to total supernova 
rate; (ii) variation of gas mass contained in a star-forming region. While mass loss from BCDs was previously 
thought to be the major cause for the variance of T>, we suggest that the other two effects are also important. We 
finally discuss the intermittent star formation history, which naturally explains the large dispersion of dust-to-gas 
ratio among BCDs. 
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1. Introduction 

Dust grains absorb stellar ultraviolet-optical light and 
emit far-infrared (FIR) light, thereby affecting the spec- 
tral energy distributions of galaxies (e.g., Takagi et al. 
|199S ) . Since the spectral energy distribution is frequently 
analysed to infer the star formation history (SFH) that 
provides us with a key to understand the evolutionary 
history of galaxies, the research on the origin of the in- 
terstellar dust is important to this issue. A key quantity 
concerning grains is the dust-to-gas mass ratio, T>. Oort & 
van de Hulst (1946) have observationally shown a strong 
correlation between the densities of gas and dust. This in- 
dicates that dust traces dense environments, which should 
be rich in heavy elements. Thus, it is worth examining the 
evolution of T> in the context of the chemical evolution of 
galaxies. 

The condensation of heavy elements is an important 
process for the formation of grains. One of the environ- 
ments where condensation takes place is the atmosphere 
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of a cool giant star (Hoyle & Wickramasinghc 1963| ). High- 
dispersion spectroscopic observation of C2 in post AGB 
stars may indicate the condensation process on C2 as a 
dust kernel (Crawford & Barlow 2000; Kameswara Rao & 
Lambert 2000). Another environment for condensation is 



a supernova (SN). Dwek & Scalo (1980) have shown that 



SNe can be the dominant source of dust grains. This is 
partly because of a rich metal content in SNe. Indeed, a 
significant amount of dust is observed within hot SN rem- 
nants (Dwek et al. 1983 ; Moseley et al. 1989; Kozasact 
al. 1989) and theoretical work by Todini & Ferrara fl200l| ) 
has explained some principal features of dust formation in 
SN 1987A. Hirashita (1999a) has shown that observational 
dust amounts suggest that ~ 10% (the fraction in mass) of 
the heavy elements ejected from stars condense into dust 
grains. Since his model prediction reproduces the observed 
trend between dust-to-gas ratio and metallicity of dwarf 
galaxies, we adopt this fraction for dust condensation. 

The processes of dust destruction should also be con- 
sidered. As shown in Dwek & Scalo (1980), dust grains are 
not only made from heavy elem ents but also destroyed in 
SN shocks (see also McKee 1989 ; Jones et al. 1996) in a cy- 
cle of the birth and death of stars. In short, T> of a galaxy 
reflects its SFH via the regulation of dust formation and 
destruction. To investigate what determines the value of 
2?, investigation of a star-forming galaxy is the most inter- 
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csting since SNe are expected to affect most largefy their 
dust amount. 

One of the observational features which give us a key 
to understand the regulation of V is the variance of T> 
itself (e.g., Lisenfeld & Ferrara 1998| , hereafter LF98). In 
this paper we examine the variance of T> among galaxies. 
According to LF98, there is a large variance in V for a 
sample of blue compact dwarf galaxies (BCDs). By defi- 
nition, galaxies categorised as BCD more or less show ac- 
tive star-forming activity. Therefore, dust formation and 
destruction are expected to occur in BCDs. Since grains 
are composed of heavy elements, application of the the- 
ory of galactic chemical evolution is useful and interesting. 
Then, we particularly consider the variance of T> in BCDs 
by using the ch emical evolution model developed by LF98 
and Hirashita ( |l999b[ hereafter H99). 

LF98 applied a chemical evolution model to explain T> 
of dwarf galaxies including BCDs. They have suggested 
that if the dust-to-gas ratio in outflow (galactic wind) is 
different from that in the interstellar medium (ISM), the 
large dispersion of T> of dwarf galaxies can be explained. In 
this scenario, a significant outflow of gas is indispensable 
to explain the observed variance of dust-to-gas ratio of 
BCDs. However, is the variance of T> determined simply 
by mass outflow? In this paper, after further investigation 
in the framework of LF98, we answer this question and 
point out that other factors are also important to explain 
the value and variance of T>. 



Indeed, a recent analysis by Tajiri & Kamaya (2002) 
has suggested that outflow is not so efficient for BCDs. 
They estimated the current momentum supply from SNe 
by using Ha luminosity, and concluded that the supplied 
momentum is not sufficient to blow away the H I envelopes 
surrounding star-forming regions and that BCDs do not 
currently suffer significant mass loss. Legrand et al. (2001) 
also suggested that low density halos around BCDs can 
be an obstacle for the ISM to escape from the galaxies 
themselves. Thus, it is worth examining mechanisms other 
than mass outflow. Since Tajiri & Kamaya ( 2002] ) adopted 
a sample in Sage et al. ( 1992 ), we also use the sample. 
Moreover, the sample in LF98 is also included because 
this paper is an extended study of LF98. 

This paper is organised as follows. First, in Sect. || we 
explain the model that describes the evolution of dust con- 
tent in a galaxy. In Sect. ||, we consider dust destruction, 
which is the most important process in this paper. Then, 
in Sect. |], model predictions are presented in compari- 
son with observations. In Sect. |^, we discuss the results 
and propose a physical mechanism that can explain the 
observations. Finally, we summarise the contents of this 
paper. 

2. Model description 

In order to consider the dust formation and destruction, 
we analyse the dust-to-gas ratio along with the chemical 
evolution model by H99. The mode l is b ased on Eales & 
Edmunds ( |l996| ), LF98 and Dwek ( |l998| ). In our model, 



we do not need to model any SFH. This has an advantage 
in considering the dust-to-gas ratio of BCDs since the SFH 
of a dwarf galaxy is generally complex (e.g., Grebel 2001 ) 
and is difficult to model. We focus especially on dust de- 
struction, because LF98 have not fully considered it. As 
we see later, our model is a powerful tool to know the 
metallicity level where dust destruction becomes effective 
enough to suppress the dust-to-gas ratio. 

2.1. Brief review of the model 



In order to investigate dust content in a galaxy, H99 has es- 
tablished a set of model equations describing dust forma- 
tion and destruction processes. In H99, a galaxy is treated 
as one zone to focus on the quantities averaged over the 
whole galaxy. The galaxy is assumed to be a closed system; 
that is, mass inflow and outflow are not considered. If the 
metallicity of the infalling material is zero or much lower 
than that of the ISM in the galaxy, the relation between 
dust-to-gas ratio and metallicity, with which we will be 
concerned in this paper, is not altered by infall (Edmunds 
[2001 ; Hirashita |2001 ). This is because the infall dilutes 
both metallicity and dust-to-gas ratio at almost the same 
rate. Our model does not include the effect of outflow, 
and this is different to LF98, in which outflow is essential 
to explain the observed variance o f the dust-to-gas ratio 
in BCDs. Since Tajiri & Kamaya (|2002|) and Legrand et 
al. (2001) have suggested that outflow is not efficient for 
BCDs, it is worth examining a case of no outflow. Indeed, 
we present another clear possibility to explain the large 
scatter of T> among BCDs later. 

The model equations in H99 (see the paper for de- 
tails; see also LF98) describe the evolution of total gas 
mass (M g ), the total mass of metals (both in gas and dust 
phases) labeled as i (Mi; i = O, C, Fe, etc.), and the mass 
of metal i in a dust phase (McLi)- We neglect dust growth 
in clouds, since Hirashita ( |1999a ) has shown that this pro- 
cess in low-metallicity systems such as dwarf galaxies is 
much less efficient than the formation of dust around stars. 
Then, we adopt an instantaneous recycling approximation 
as in L F98 and H99 according to the formalism in Tinsley 
( |1980D : Stars less massive than m t (present turn-off mass 
set to be 1 M Q ) live forever and the others die instanta- 
neously. 



2.2. Solution of the model 

Dust-to-gas ratio and metallicity of galaxies are observa- 
tionally known to correlate with each other (e.g., Issa et 
al. 199C). This relation has recently been used as a test for 



chemical evolution models including dust formation and 



destruction (LF98; H99; Hirashita p99a| ; Edmunds |200lD 
The model by H99 reduces the following differential equa- 
tion: 



y, 



(1) 
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where T>i is the mass fraction of an element i locked up in 
dust (i.e., T>i = Md,i/M g ); Xi is the mass fraction of an 
element i (i.e., Xi = Mj/Mg); /i n ,i quantifies what frac- 
tion of an element i ejected from stars condenses into dust 
grains; /?sn is "dust destruction efficiency," which we ex- 
plain later; 1Z is the fraction of stellar mass subsequently 
returned to the interstellar space, and 3-j is the mass frac- 
tion of an element i newly produced and ejected by stars.^ 
The definition of /?sn is as follows: 



0s* 



tsn-0 



(2) 



where tsn is the timescale of dust destruction by SN 
shocks (Eq. |5|), and ip is the star formation rate (SFR). 
Hereafter, /?g N is called "dust destruction efficiency", be- 
cause it is inversely proportional to the timescale of dust 
destruction (tsn)- This is the most important parameter 
in this paper. 

When all the quantities except Xi and T>i are constant 
in time, the analytical solution obtained by LF98 is appli- 
cable. With our notations, it is rewritten as 



Vi(Xi) = -X t + (1 
a 



(3) 



where a = (11 + ^sn)/^, b = f m .iU/yi, and c = f in .i. 

Here, we select oxygen as a traced element (i.e., i = O) 
according to LF98, because (i) most of the oxygen is 
produced by massive stars (Type II SNe and their pro- 
genitors), (ii) oxygen is one of the main constituents of 
dust grains, and (iii) the common tracer for the metal 
abundance in BCDs is an oxygen emission line. The first 
item (i) means that an instantaneous recycling approxi- 
mation may be reasonable for the investigation of oxygen 
abundances, since the generation of oxygen is a massive- 
star- weighted phenomenon. In other words, results are 
insensitive to the value of m t . Following H99, we adopt 
(K, y ) = (0.32, 7.2 x 1(T 3 ), which are consistent with 
the relation between dust-to-gas ratio and metallicity of 
nearby galaxies including our BCD sample. 

3. Destruction efficiency of dust 

In this paper, we discuss the variance of dust-to-gas ratio 
of a BCD sample in terms of the variation of /3 sn de fined 
in Eq. (||). We estimate /3sn based on McKee (198S) and 
LF98, while we also address their differences. 

We assume that gas is divided into two components: 
gas in the star-forming region and that in the H I en- 
velope. Such an envelope is generally observed around a 
star-forming region of a BCD (e.g., van Zee et al. 1998[ ). 
We denote the gas mass fraction in the star-forming re- 
gion as Xsf and that in the H I envelope as Xm (i-e., 
Xsf + Xhi = 1). We distinguish the two regions for the 
comparison with the IRAS sample (Sect. O). IRAS FIR 



1 TZ = R and y t = y(l - R) for the notation in LF98. These 
two parameters are calculated in the formalism of the instan- 
taneous recycling approximation. 



bands are sensitive to dust hotter than about 25 K. Such 
"warm" dust exists in star- forming regions, not in H I en- 
velopes. Calzctti et al. (1995) have also shown by using the 
IRAS sample of actively star-forming galaxies that 70% of 
the FIR flux comes from such a warm component of dust. 
Thus, dust mass derived from the IRAS observation of a 
BCD is considered to trace the dust in the star-forming 
region. This suggests that it is useful for us to consider 
dust contained in star-forming regions as long as we are 
interested in the comparison of our result with the IRAS 
observations. 

Thus, we estimate tsn and /?sn in star-forming regions. 
Gas mass accelerated to a velocity of v s by a SN, M s (v s ), 
is estimated as 



E. 



M s {v s ) = 6800^- M c 



(4) 



where £51 is energy released by a SN in units of 10 51 erg 
and v S 7 is v s in units of 10 7 cm s _1 (McKee 1989). In 
swept ISM, dust is not fully destroyed. Thus, the fraction 
of destroyed dust, e (~ 0.1; McKee 1989), should be mul- 
tiplied. Since we are interested in mass swept by multiple 
SNe in a star-forming region, SN rate must be considered, 
which is denoted as 7. Then, tsn is expressed as 



M g X SF 



e 7 M s (100 km s" 1 ) 



(5) 



Here, we substitute v s with the threshold velocity for dust 
destruction (100 km s _1 ; McKee 1989). Comparing equa- 
tions (j|) and (||), we obtain 



/3 S n = eM s (100 km s" 1 ) 



1 



7 



(6) 



We should consider both Type la and II SNe in esti- 
mating 7. Type II SNe occur "soon" after a star forma- 
tion, because the progenitors of Type II SNe have lifetimes 
much shorter than the age of the universe. On the other 
hand, Type la SNe, whose progenitors are low-mass stars, 
occur about 1 Gyr after a star formation. Thus, 7 and /?sn 
should be sensitive to the SFH. We express 7 as 



7 = 7la + 711 



(7) 



where 7i a and 711 are the rate of Type la and Type II 
SNe, respectively. We note that LF98 estimated /?sn to 
be 5. But, for example, the fraction between 7 i a and 711 
changes according to SFH of BCDs. In the next section, we 
examine the relation between dust-to-gas ratio and metal- 
licity for various /3sn, whose probable range is constrained. 

4. Sample and result 

4.1. Observational samples 

We compare the model calculation in Eq. (||) with the ob- 
served values of dust-to-gas ratio and metallicity of BCDs. 
We a dopt th e data in Table 2 of LF98 and Table 1 of Sage 
et al. (1992). The latter sample has been adopted by Tajiri 
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Table 1. Data of Blue Compact Dwarf Sample 



Object D S eo S 10 o log M m log M' d RAS 12 + (O/H) ref. a 

(Mpc) (Jy) (Jy) [M ] [M ] 



II Zw 40 


9.1 


6.6 


5.8 


8.30 


4.55 


8.15 


S92 


Haro 2 


20.3 


4.8 


5.5 


8.68 


5.39 


8.4 


S92 


Haro 3 


13.7 


5.2 


6.7 


8.76 


5.21 


8.3 


S92 


UM 439 


12.6 


0.39 


1.2 


8.23 


4.96 


7.98 


S92 


UM 462 


11.9 


0.99 


1.1 


8.15 


4.21 


7.89 


S92 


UM 465 


13.2 


0.99 


1.3 


7.71 


4.48 


8.9 


S92 


UM 533 


10.4 


0.51 


0.54 


7.76 


3.75 


8.10 


S92 


UM 448 


6.00 


4.14 


4.32 


7.54 


4.16 


8.08 


LF98 


IC 3258 


21.2 


0.490 


0.970 


8.55 


5.03 


8.44 


LF98 


Mrk 7 


42.3 


0.480 


0.970 


9.56 


5.64 


8.54 


LF98 


Mrk 33 


21.6 


4.68 


5.30 


8.77 


5.41 


8.40 


LF98 


Mrk 35 


14.5 


4.95 


6.74 


8.73 


5.29 


8.30 


LF98 


Mrk 450 


12.1 


0.480 


0.820 


7.77 


4.37 


8.21 


LF98 


NGC 4670 


12.1 


2.63 


4.47 


8.52 


5.10 


8.30 


LF98 


NGC 4861 


12.9 


1.97 


2.26 


9.13 


4.60 


8.08 


LF98 


II Zw 70 


17.6 


0.710 


1.24 


8.56 


4.89 


8.11 


LF98 



a S92 — Sage et al. (1992); LF98 — Lisenfeld & Ferrara (1998) 



& Kamaya ( [20021) . We select BCDs whose 21 -cm H I emis- 
sion and FIR dust emission are both detected. The data 
are summarised in Table [j], 

The observational dust-to-gas ratio in a star-forming 
region should be 



V 



obs 



Mi RAS /(M m X SF ) , 



(8) 



where M^ RAS and Afjji are dust mass determined from 
the IRAS observation and the total H I gas determined 
from the 21-cm observations. Mm traces H I gas in both 
the star-forming region and the envelope of a galaxy. On 
the other hand, M^ RAS traces the amount of dust in the 
star-forming region, because the IRAS is sensitive to high- 
temperature (il 25 K) dust. In the following two para- 
graphs, we describe the two quantities further. 

LF98 divided the observed H I mass by a factor of 2 
to obtain the gas mass in only the star-forming region. 
In other words, LF98 assumed that X SF is 0.5. However, 
Xsf is hardly constrained observationally. In this paper 
we assume Xsf = 1 for the observational sample to obtain 
the first result. As we discuss in Sect. |5|, the variation of 
Xsf also contributes to the variation of /3sn- Therefore, 
Xsf affects our analysis in the following two ways: 

1. Observationally, the estimate of dust-to-gas ratio in a 
star- forming region is affected by the value of Xsf- 

2. Theoretically, the value of Xsf affects the value of /3sn, 
thus changing the dust-to-gas ratio. 

It is difficult to discuss the former quantitatively, because 
we need a high-resolution spatial map of H I distribution 
and a reasonable observational definition of a star-forming 
region on the map. However, we can discuss the latter 
because the effect of Xsf on /3sn is well determined from 
Eq. #. 



The observational dust mass, M^ RAS , is derived from 
the luminosity densities at wavelengths of 60 fim and 100 
/^m observed by IRAS using Eq. (4) of LF98. The IRAS 
bands are insensitive to the cold (Ss 20 K) dust that lies 
out of star-forming regions. Moreover, since the dust in 
a star-forming region suffers destruction by SN shocks, 
we should take into account an efficient destruction in 
the star-forming region. Therefore, in order to discuss the 
IRAS sample and the selective dust destruction in star- 
forming regions, we need to define /?sn as the dust de- 
struction efficiency (Sect. |3|) in star- forming regions as we 
have done in Eq. (||). 

Finally, we hypothesise that the fraction of oxygen con- 
tained in dust grains is constant for all the sample BCDs. 
Following H99, we assume the Galactic composition of the 
grains: 



V = 2.2V, 



O • 



(9) 



4.2. Conclusion from our model 



In Figure [j} we show analytical results calculated accord- 
ing to Eqs. (||) and (||) with i = O for various values of 
Psn as solid line (/?sn = 1), dotted line (/3sn = 5; the case 
of LF98) and dashed line (/?sn = 25), respectively. The 
black and gray squares indicate the observational samples 
in Sage et al. ( |1992| ) and LF98, respectively. The num- 
ber ratio of oxygen atoms to the hydrogen atoms is de- 
noted as (O/H). We convert the mass fraction of oxy- 
gen, Xo, to (O/H) for the model prediction, assuming 
log X + 10.80 = 12 + log (O/H). In this figure, / inj0 = 0.1 
is assumed to concentrate on the variation in /3sn- LF98 
have shown a large variety of f lni q to explain the large 
scatter of the relation between V and Xq in their Figure 
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Fig. 1. Relation between dust-to-gas ratio and metallicity. 
The solid, dotted, and dashed lines show the model results 
with the dust destruction efficiencies (Eq. |2|) /3sn = 1, 5, 
and 25, respectively. The black and gray squares show the 
data for the blue compact dwarf galaxies in Sage et al. 
(|1992|) and in Lisenfeld & Ferrara ( 1998 ), respectively. 



7. Since Hirashita (1999a) has shown that /; n . o ~ 0.1 
reproduces the observed trend between dust-to-gas ratio 
and metallicity, it is worth considering the effect of /3sn 
by setting / in ,o = 0.1. 

From Figure |l|, we see that the variance in dust-to-gas 
ratio is reproduced by an order-of-magnitude variation in 
/3 SN (1-25 here). We note that LF98's value (/3 S n = 5) is 
within this range. Even if outflow does not efficiently occur 
in BCDs, we can explain the variance of the dust-to-gas 
ratio of BCDs with the various "destruction efficiency" of 
dust, /3sn 7 although we can never reject the importance 
of outflow in the framework of this paper. The important 
point is that we have demonstrated that dust destruction 
by SNe can play an important role in producing the vari- 
ance of V. The result also indicates that the difference 
in dust-to-gas ratio among the three lines becomes clear 
at an oxygen abundance, 12 + log(0/H) ~ 8. Thus, we 
conclude 

1. that the dust destruction by SNe can vary the dust-to- 
gas ratio when the metallicity, 12 + log(0/H), reaches 
about 8 (~ 10% of the solar metallicity), and 

2. that the large variation of the dust-to-gas ratio among 
BCDs is explained by the variation of dust destruction 
efficiency, /?sn- 



5. Discussion 

5.1. Scenario 

We propose a scenario for large variation of /?sn for each 
BCD along with intermittent SFH. The time variability 
of SFH on a short (< 1 Gyr) timescale has been sug- 
gested observationally by Searle & Sargent (1972) and 



theoretically by Gerola et al. (1980). Nonlinear processes 
in the ISM may also cause an intermittent star formation 
(Ikeuchi 1988; Kamaya & Takcuchi 1997 and references 
therein). Thus, it is pertinent to consider the intermittent 
SFH of BCDs. 

As shown in the following, an intermittent star for- 
mation history leads to the time variation of /?sn, because 
j/ip depends on time. We consider an intermittent SFH: a 
starburst whose SFR is "0burst and an inter-starburst epoch 

whose SFR is Winter- We aSSUme that ^burst = lOO'i/'inter, 

for example. Such a two-orders-of magnitude variation in 
SFR is proposed theoretically by Gerola et al. (1980) and 
Kamaya & Takeuchi (1997). While the starburst is go- 
ing on, we expect that almost all the SNe are Type II 
(Tburst ~ Til, where 7burst is a typical SN rate in the 
bursting epoch). However, in the inter-burst epoch, Type 
la SNe can be dominated (7i n ter ~ 7ia, where 7i n tor is 
a typical SN rate in the int er-bur st). According to the 
model by Bradamante et al. ( 1998 ), a given stellar popu- 
lation releases energy in the form of Typ e II a nd la SNe 
with a ratio of 5:1 (Bradamante et al. 1998 ; see their 
Fig. 9. This value is essentially determined by the ini- 
tial mass function (IMF), and they assumed the Salpeter 
IMF with a stellar mass range from 0.1 to 100 M Q f\ 
Since they assumed the same energy between Type la 
and II SNe, this means that the number ratio between 
Type la and II SNe is 5:1. Therefore, we expect that 
7inter ~ 7burst/5. The intermittent star formation finally 
predicts that 7 inter /Winter ~ 207 bu rst/V'bur S t- This means 
that an intermittent SFH can cause a 20-times variation 
in /?sn oc j/ip during a single star formation cycle. 

Furthermore, Fig. [l] shows that the value of /3sn has 
little effect on the relation between dust-to-gas ratio and 
metallicity for 12 + log (O/H) < 8. Therefore, until the 
metallicity level becomes 12 + log (O/H) ~ 8, the relation 
between dust-to-gas ratio and metallicity evolves in the 
same way whatever the value of /3sn might be. On the 
contrary, the relation is largely affected by /?sn if 12 + 
log (O/H) > 8. Then, we study the response of the relation 
between dust-to-gas ratio and metallicity to the change of 
/3sn at 12 + log (O/H) ~ 8, as we are interested in the 
intermittent SFH. 

First, we shall estimate a typical metallicity increment 
during a single star formation epoch of the intermittent 
SFH. The metallicity increase during an episode of star 
formation, AZ, can be estimated by AZ ~ yM*/M s , 
where A/* is the mass of stars formed in the episode, and 
y is a chemical yield. If the IMF is similar to that of the 
Galaxy, y ~ Z Q (i.e., AZ -> Z Q for M* -> M g ). We es- 
timate M» by multiplying observed SFR with a duration 
of an episode of a star formation activity. Assuming that 
the SFR is 0.1 _M© yr -1 and that the duration is 10 7 yr 
(Legrand et al. |200l| ), we obtain A/* ~ 10 6 M©. With typ- 
ical gas mass M g ~ 10 7 Mq, we obtain AZ ~ O.IZq. This 



2 We assumed the Salpeter IMF with the range from 0.1- 
120 Mq, but the difference in the higher mass cut has little 
influence on the SN rate. 
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Fig. 2. Relation between dust-to-gas ratio and metallicity. 
The three dotted lines show the model results same as Fig. 
|l|. The squares are the data points same as Fig. |l|. The two 
solid lines represents the result of the calculations which 
change /3sn (from 5 to 1 and from 5 to 25 for upper and 
lower lines, respectively). 



corresponds, for example, to the metallicity increase from 
12 - flog ( O/H) = 8.0 to 8.2. The model by Bradamante et 
al. ( 1998| ) also indicates that one episode of star formation 
can result in such a metallicity increment. 

As shown above, the effect of intermittence can be ex- 
amined by changing /?sn- In order to examine the effect 
of time variation of /3snj thus, we calculate the relation 
between dust-to-gas ratio and metallicity in the following 
two cases: 



1. /3 SN = 5 for 12 4 
12 + log (O/H) > 

2. /3 SN = 5 for 12 - 
12 + log (O/H) > 



log (O/H) < 
log (O/H) < 



and /3sn = 25 for 
i and /3sn = 1 for 



In Figure ||, we show the result of the two calculations (two 
solid lines). The lower and upper branches represent the 
cases 1 and 2, respectively. The three dotted lines show 
the results as in Fig. [l| We see that when the metallicity 
12 4- log (O/H) increases from 8.0 to 8.2, the line rapidly 
converges to the two dotted lines, which represent the re- 
sult for constant /?sn (1 and 25, respectively). This rapid 
convergence supports the idea that the dust-to-gas ratio 
varies widely in response to the time-evolution of /3sn in 
an episode of star formation at the metallicity level of 
12 4- log (O/H) ~ 8. Thus, it is possible to determine the 
variance of the dust-to-gas ratio among BCDs from the 
time variation of /?sn- 

5.2. Age difference 

If the age of BCDs varies, the present turn-off mass of 
stars is changed. As a result, the returned fraction of gas 
(TZ), the metal yield (D^), and the dust supply from stars 



(/in,i) arc effectively different among BCDs. The effect of 
varying 1Z and ^ on the relation between dust-to-gas ratio 
and metallicity has been examined by H99.[] However, the 
resulting relation is less sensitive to the two parameters 
than to /3s n- The dependence of fin,i on the turn-off mass 
can be important since it largely affects the dust amount 
in low-metallicity systems (LF98; Hirashita 1999a). 

Thus, if the BCD sample proves to have a large age 
variation, we should reconsider the variation in dust-to- 
gas ratio with a time-dependent formulation. Indeed, we 
cannot reject the possibility that some BCDs are much 
younger than the cosmic age. A metal-poor BCD SBS 
0335- 052 may be younger than 5 x 10 6 yr ( Vanzi et al. 
fcOOOp . Recently, Hirashita, Hunt, & Ferrara (f(X)§) have 
succeeded in explaining the dust amount of SBS 0335-052 
with a time-dependent formulation applicable to young 
galaxies. 

5.3. Future observations 

The time variation of Agp also leads to the time depen- 
dence of /3gp (Eq. ||) , although Xgp was assumed as unity 
(Sect. 4.1). If the various Xgp for the BCD sample is inter- 
preted to reflect the time evolution of Xsf in each BCD, 
we can suggest that the gas mass in a star-forming region 
should change temporally because of the mass exchange 
between the star-forming region and the envelope. Such 
a mass exchange during episodic star formation a ctivit y 
in BCDs is indeed suggested by e.g., Saito et al. ( 2000 ). 
The temporal change of Asf is also possible if the ISM 
in the star-forming region is consumed for star formation 
and locked in stellar remnants like white dwarfs, neutron 
stars, and black holes. 

In order to constrain ATsf , we need to observe H I emis- 
sion or FIR emission with angular resolution fine enough. 
The present typical angular resolution of 1' corresponds 
to 2.9 kpc in physical size if a galaxy lies at a typical 
distance of 10 Mpc. Future large space FIR telescopes 
such as Hersche m (e .g., Pilbratt |2000| ) or SPICA fj (e.g., 
Nakagawa et al. 2000D will resolve the star-forming regions 
of the BCDs. For example, the Japanese future infrared 
satellite SPICA will have a diameter larger than 3.5 m. If 
the diffraction limit is achieved, the angular resolution be- 
comes 6" at the wavelength of 100 /im. This corresponds 
to 290 pc at the distance of 10 Mpc, and is comparable 
to or smaller t han th e half-light radius of a typical BCD 
(Marlowe et al. |l999 ). 



6. Summary 

In this paper, we have analysed the relation between dust- 
to-gas ratio and metallicity of BCDs by using our chemical 

3 See Fig. 3 of the paper. Although H99 examined "differ- 
ent IMF", this difference only affects the values of 1Z and 3^- 
Therefore, his result can be used to see the effect of varying TZ 
and Vi. 

4 http:/ /astro. estec.esa.nl/First/ 

5 http:/ /www. ir.isas.ac.jp/SPICA/index. html 
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evolution model. We have focused on the dust destruction 
process, because this process was not investigated in LF98. 
We have shown that the dust destruction significantly af- 
fects the dust-to-gas ratio when the metallicity is larger 
than 12 + log(0/H) ~ 8. The intermittent SFH can ex- 
plain the large variety of dust-to-gas ratio among BCDs 
through the "dust destruction efficiency" of /3s n- 
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